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PROBLEM 

In 2009, the U.S Strategic Command space objects (assets and debris) tracking update: 

• ~ 19,000 space objects (diameter >10 cm) are being tracked 

• > 600,000 space objects (10cm >diameter>l cm ) are unobserved 

Recent impacts have also added more objects (space debris) to track: 

• Satellite Collision Russian Cosmos 2251 and the Iridium satellite in 2009 

• Destruction of Chinese satellite in 2007 

REQUIREMENTS 

Space objects require: 

• Collision avoidance mitigations 

• Orbit maintenance/maneuvers (assets) 

• Cataloging and identification (debris) 

Specifics for the requirements: 

• Performed within required accuracies 

• Cost effective : Applicable to large numbers of space objects 


Artist rendering of space object collision, courtesy CNES 
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CURRENTLY 

Statistical Orbit Determination approach 


Estimated trajectory 
Reference trajectory 
True trajectory 


Dynamjc model 
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X=—^r +perturby»£ accelerations 



Apriori Gaussian 
uncertainty 


x= 


X 


K 

Y 


V y 

Z 


v z 

K 

+ v(0 = 

-juX/r 3 

K 


- HY/r 3 

K 


1 

cn 

'nT 

i 

i 


v(0 



Non-Gaussiair 

uncertainty 


Prediction Period 


Measurement model 



P = yl(r -R)*(r-R) + w(t) 
p = [(r-R)»(r-R)]/p+ w(t) 
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State Update 


Ground Station 
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PROBLEM STATEMENT 

1. Need Full Probability Density Function (PDF) representation for low probability 

events present in the heavy tails (Non-Gaussian) 

2. Use a nonlinear filter that is capable of full non-Gaussian PDF state estimation 

3. Need a compressed representation of this PDF distribution for real case scenario 
applications 
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• Extended Kalman Filter (EKF) (Tapley, Shutz and Born 2004) 

- Nonlinear models assuming Gaussian uncertainties 

- Linearization about a current mean (computes Jacobian) 


• Unscented Kalman Filter (UKF) (Julier and Uhlmann 1997 et al) 

- Uses a series of weighted sample points to approximate the mean and covariance 


Actual (Sampling) Linearized (EKF) ITT^-sigma points 



http://ars.els-cdn.corn/content/image/l-s2.0-S0021999108000132-gr4.jpg 


• Batch Weighted Least Squares (Gauss 1809) 

— Estimates the state x by minimizing the performance index 




£ r We 


8 




Gaussian Mixtures Model (GMM)(DeMarsetal2011) 

l Covariance P 

p(x) = Y J w i p g (x; J u i ,P i ) 


i = 1 


0 


Detects Nonlinearities 
It splits 





Nonlinear 

Propagation \J 

Figure Ref: J .T. Hartwood et al. Gaussian Sum filters for Space Surveillance: Theory and 
Simulations . JGCD, Vol. 34. No. 6. Nov-Dec 2011 



where p is a Gaussian PDF 

- Differential entropy d/dt[H(x)] is used as a measure of nonlinear detection 

- Uses UKF as a filter for each p g Entropy 

Polynomial Chaos Expansion (PCE) (Jones et al 2012) 

- Solution of the Stochastic Differential Equations as linear expansions of multivariate 
polynomials 

- PCE estimates the coefficients c a of the expansion 



My approach: 

Propose using the Particle Filter 

• Fully captures PDF (Probability Density Function) distribution 

• Uses a large number of samples to approach the optimal estimate 
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The Particle Filter (PF) is a sequential nonlinear estimator that 

• Uses random independent particles to represent the random state x as a PDF 

• As N increases, we reach the optimal PDF representation 



• actual state value actual state trajectory 

x measured state value estimated state trajectory 

• state particle value particle propagation 

state pdf (belief) particle weight 


Figure Ref : Bhaskar Saha, Introduction to Particle Filter. Ames Research Center. April 16, 2008 
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initialization 


Time Update 


Measurement Update 


Effective # of Samples 


Sampling Degeneracy 



Resample 







Sample Degeneracy when data outliers occur or when measurement noise is small 



CIRCUMVENTION: 


• Resampling of particles (N effective samples < N threshold ) 

• Adaptive PF: Soto (2005), Gang and Xiao-Jun (2008), Hwang and Speyer (2011). 

• Other PF Methods studied: Auxiliary Particle Filter, Regularized Particle Filter, Rao- 
Blackwell/Marginalized Particle Filter 


Our Approach for PF in Statistical Orbit Determination: 

Used Generic Particle Filter 

• Effective: sporadic measurements 

• "Jitter" resampled particles by adding noise 
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'' Multivariate Density 

i\j?(X,Y,Z,V x ,V r ,V z ) 


(A) Particle Filter: 

• Meas. Updates and Predictions to B 

• Large number of particles - compression 


(B) 


Predicted State 

LxN 

(states EC x particles) 


NFA +. 


FastICA 


Uncorrelated (UC) 

LxN 

(states uc x particles) 


Binning 

> 


Probability Mass Fen 
(PI VIF) 

L x #bins X 2 
(states uc x probabilities) 


WT/FFT + 


Information 

Measures 


Compressions 


As Lossless 
as possible 


Coefficients for PMF 
LxC 

(states uc x coefficients) 


Inverse Transforms 


Reconstructed 

State 

LxN 

(states ECrec x particles) 


Inverse 


FastICA 

+..NFA 


Compressed & 
Uncorrelated (UC) 

LxN 

(states uc x particles) 


WT/FFT 


Discrete 


Sampling 


Compressed PMF 

L x #bins + 3L 
(states uc x probabilities 
x abscissa values) 
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Demonstrated the use of the PF as a statistical OD method for non-linear 
state estimation and uncertainty predictions as a full state PDF 

Used nonlinear multivariate decorrelation methods to decorrelate the 
multivariate state PDF 

Performed univariate uncorrelated state PDF compressions for data 
allocation and transmission cost reductions 

Demonstrated the potential for cost effective ephemeris storage for the 
cataloging of space objects and debris 
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The Magnetospheric Multi-scale mission (MMS) Phase I 


Normal 



Orbital elements 

a 

42095 

e 

0.8182 

i 

28.5° 

Q 

357.857° 

U) 

298.225° 

M 0 

180° 


The VBN frame is the reference path orbital velocity, binormal and normal directions. 

• Rotating frame 

• Easy to visualize propagations 
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Skewness 



Illustrating the non-Gaussian content of the PDF propagations of the MMS orbit for 1 period 
• Apogee as the initial condition 



Skewness: a measure of asymmetry 
(3 rd order moment) 


Skewness = 


[£{( X -/ y ) 2 }] 3/2 


Excess-Kurtosis: a measure of peak intensity 
(4 th order moment) 


Excess - Kurtosis 


e{(x-/j) 4 } 

[E{{X-nf}f 
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Epoch: 04-NOV-2011 16:00:00 TDRS-l and TDRS-3 



00:00:00-01:15:00 


Estimation of state vector X: using the Extended Kalman Filter(EKF). 

using the Particle Filter for 1000 particles . 


CASE 1: Over 1 orbital period with 4 measurements 

CASE 2: CASE 1 AND an additional 4 orbital periods without measurements 


PF code was developed for the Orbit Determination Toolbox (ODTBX) that uses some of its plotting 
capabilities. 21 




Transformed Position States from XYZ ECI to VNB 
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EKF PF 



Meas. update errors 
Ensemble mean error 
2 a from P formal (EKF) 
2o from P empirical (PF) 
Ensemble errors 
1-3 a of ensbl. errors 


NOTE : Resampling 
threshold N thr = 0.75 N 

PF outperforms EKF 
based on mean errors 

PF has larger 
uncertainty bounds vs 
EKF 


Transformed Position States from XYZ ECI to VNB 
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5000 
4000 
3000 
2000 
1000 
V(m)0 
-1000 
-2000 
-3000 
-4000 

0 1 2 3 4 5 


EKF 



PF 


- Ensemble mean error 
+ 2 a from P formal (EKF) 

+ 2o from P empirical (PF) 
| Ensemble errors 

| 1-3 a of ensbl. errors 


Max. Prediction errors 
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200 | 1 . 1 1 1 


Post-updates 0.25N thr 



PF uncertainty 
bounds are tighter 
than EKF 

Bi-normal component 
has good accuracies 
at both apogee and 
perigee for the PF vs 
EKF. 


Transformed Position States from XYZ EO to VNB 1 
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-2°° o J ^ ^ -J J .200 1 l_ 

Transformed Position States from XYZ EC }to VN& 


- Ensemble mean error 
+ 2 a from P formal (EKF) 

+ 2o from P empirical (PF) 
| Ensemble errors 

| 1-3 a of ensbl. errors 


Max. Prediction errors 



• Post-updates 0.75N thr 

• PF uncertainty 
bounds are tighter 
than EKF 

• Bi-normal component 
has good accuracies 
at both apogee and 
perigee for the PF vs 
EKF. 
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Epoch: 04-Nov-2011 16:00:00 



Direction 


Particle Filter for 1000 particles vs. GMM using UKF 

Scenario: Over 1 orbital period with 15 min measurement update at end of 1 period 

View the State PDF at the Measurement Update instance: A) As planar contour plots 

B) As 3-D PDF plots 
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z-Position y-Position (km) 



71768.842 


GMM + UKF 


1000 MC data points 
^^Probability Contours 

♦ GM 1 (0.2236) 

♦ GM 2 (0.5367) 

♦ GM 3 (0.2397) 
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6963.3 



1 1 


1000 Particles 
Cfi>Probability Contours. 
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GMM + UKF 


x 10 


+ GM 1 (0.2236) 
GM 1 samples 
GM 2 (0.5367 
GM 2 samples 
♦ GM 3 (0.2397) 
GM 3 samples 



0 

71768.850 
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71768.820 

y -Position (km) 


x 10" 


6963.3 


6963.35 6963.4 

x-Position (km) 


6963.45 




#GM1 ( 
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— $GM 2 ( 
•® GM2 s 
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— e GM3 ! 
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PF Probability Scale 
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0 
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: p<^\0.0006 
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PRINCIPAL COMPONENT ANALYSIS (PCA) 

x : zero mean data 

V =WAW t . . 

^ x ” 1 '- ,r s: principal components 

s = W T x 


PCI 




i 

PC 2 

1 

>C ^ 

.... . , 


PCA finds the principal components in a new transformed coordinate system with the greatest variance as 
the first component and the subsequent ones follow. 

INDEPENDENT COMPONENT ANALYSIS MCA) 

X = Bs x : whitened data 

s : independent states 

s-B t x 

ICA finds the linear combinations of the unit 
sphered observations x (w T x) that are in the 

direction w of the max or min Kurtosis (4 th order cumulant) B = {wl w2}. 



Unit sph ere/ circle 


FastICA : uses a fixed point iteration scheme (MUCH FASTER). 

Kurtosis^.) = is {s/} -3(E{sf}) 2 where ^ = w T x 

w(i + 1) = E {x(w(i) T xf } - 3w(i) iterated |w(/ + l) r w(/)|— >1 
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Orbit parameters: eccentricity = 0.2, P 0 ^ o 2 w = (1000m) 2 and o 2 vx vyvz = (lm/s) 2 

1 rev of propagation, Epoch of interest at perigee. Period = 18hrs 


• Initial state dimension L = 6 icaotpca 4 components 

• Choice of 4 : Out-of-plane motion is decoupled from in-plane motion 

- Fundamental nonlinearity comes from Kepler's equation 

• Demonstrate potential for augmented states not necessarily required for accurate state 
prediction (Station location errors, Range biases etc) 


• The reconstructed distributions are measured using: 

• The Kolmogorov-Smirnov (K-S) test 

— Quantifies the max. distance D between 

the two Cumulative Distribution Functions (CDF) 

D ab = sup \F a (X) -G b (X)| 

X 



• The state particles are scaled to canonical units 


http://www.statsoft.eom/Portals/0/Support/ks%20graph.JPG 


• Ensures equal weightings for velocity components (affects eigenvalues) 

• Distance : (1DU) = 6378.145km, Velocity : (1DU/ TU) = 7.90536828km/s 


31 
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POSITION: X-Y plane (DU) 


ICA 
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□ Original Data 
• Reconstructed Data 


Excess 

Kurtosis 

PCA 

ICA 

1 st 

1.30 

24.39 

2nd 

22.66 

2.55 

3 rd 

0.46 

0.02 

4th 

-0.08 

-0.4 


K-S Statistic 

PCA 

ICA 

X 

0.0740 

0.0080 

Y 

0.0100 

0.0020 

Z 

0.0030 

0.0150 


Shortcomings of PCA/ICA: 

* For high non-Gaussian 
behavior the components 
were not fully 
decorrelated 

Solution 

* Use Nonlinear Factor 

Analysis ^ 






The Nonlinear Factor Analysis (NFA) is a nonlinear mapping of the sources s to the 
observations x modeled by a multilayer perceptron (MLP) network model. 


x = f(s,0) + n 

= %(j){/4s + a) + b + n 

= % tanh(^ + a) + b + n 



Our goal is to estimate the sources s and the unknown variables 0 =(/^,^,a,b) 

A E y{ H n* L , % E y{ LxHn , a e y\ !/ " and b E ${ L based on minimizing the cost function 
between the posterior p and the approximate q given by the Kullback-Leibler Divergence 
metric D 

D(q(& \c)\\p(&\x,'»))= | q(& 1 £) In d® wtlere 

p(® | x, »?) are ran £ om varia6ks 

^ and » are some prior assumptions. „ 


Ref : Theory and MATLAB code developed by Dr. Annti Honkela at Helsinki University of Technology, Finland 




The prior assumptions ^ and # are given as follows: 


x - S tanh( / ^s' + a) + b + n 
x ~ N(f(s, 6), diag(e 2r " )) 
s ~ N(0,diag(e 2V ’)) 

Parameters 

A, B,a,b~ N(m e , diag(e 2Vd )) 

0 = {A,B,d,b,I} and 

0 = {^ 4 , 5 , 5 , 5 ,£} 

are the solve-for parameters. 

Known: x and priors for 0 and © 
and 

Hyperparameters: 


The variance is parameterized to v = log a 
D = f q(& | g )ln g(Q ^ g) f/0 

/>(© | X, *?) 

* f(0) 

= fg(0|^)ln J0 + lnj^f^^ 

/>(0, x | W ^ 

C = fg(0|£)ln </© 

/?(0, x | fW 

C = ^[lng(0|^)] + ^[-ln/7(0,x|^)] 
min C = C q + C p = E[ In g(0 1 £)] + £[- In /?(0, x | #)] 
w.r.fc 0 and 0 


^£ = S + S=o 


dc dc dc 


_ q 


+ ■ 


p _ 


= 0 


v„,v s ,v Bj ~ N(m v ,diag(e 2K )) 

(m a , v a , m b , v b , m v _,v V ',m VB ,v Vii ) ~ MO, 100" (uninformative priors 


50 50 50 50 50 50 

xv Here & and 0 t/ie mom and variance respectiveCy 

2 \ • r .. • 34 




The obtained mean and variance terms for s and 0: 

• mean s: L x N 

• variance s: L x N 

• mean 0:^(H n x L), H n ), a(H n x 1) and b(Lx 1) 

• variance 0: ^(H n x L), x H n ), a(H n x 1 ) and b(L x 1 ) 


L : state dimension 
N: No. particles 
H n :No. Neurons 


The variance terms for the sources s are assumed to be zero (and variances 0=0) 



• FastICA produces uncorrelated sources s' 


s'l s'N 


> 
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The obtained mean and variance terms for s and 0: 

• mean s: L x N 

• variance s: L x N 

• mean 0:^(H n x L), H n ), a(H n x 1) and b(Lx 1) 

• variance 0: ^(H n x L), x H n ), a(H n x 1 ) and b(L x 1 ) 


L : state dimension 
N: No. particles 
H n :No. Neurons 


The variance terms for the sources s are assumed to be zero (and variances 0=0) 



• FastICA produces uncorrelated sources s' 


s'l s'N 


> 
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The uncorrelated sources s', are binned into a 2 7 histogram then transformed into a 
normalized discrete probability density function (PMF) p for compression. 


FAST-FOURIER TRANSFORM 

• The FFT calculates the coefficients of the discrete PDF as follows: 

A — t/-ix*-i) 

c(k ) = 2^p(J)e n where (k = 1, n) p(j) is the discrete probability value 

j = i normalized from the histogram 


WAVELET TRANSFORM 


• Use a wavelet filter (Daubechies 2) to extract terms that represent the PMF p as a 

function of approximate and detail coefficients ^ 

«,<,["] = 2>M*M2h-A:] 

k 


P(k) 



7 ,[n] = ^p[k]* g i2n-k] 

k 

00 

pW=Z((“yo[»]*M2n-^])-Kr ; [n]*g[2»-^])) 

n =— oo 

y'O is the first level of decomposition 

j are the subsequent levels of decomposition 
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The uncorrelated compressed PMF representations are reconstructed by taking the 
inverse transforms of the FFT and the WT. 


I FFT: p{j) = ^c{k)e n where (7 = 1 , n) and (m<n) 

k=\ 

F 

IWT: p(k)= 'Y' J ((oc i g[n]*h[—k + 2n])+(y i [n]*^-k + 2n]y) where (k = 1, n) 

n=-F 



• Filters (WT) 

• Abscissa values 

• Coefficients 


Discrete sampling . Mixing matrix A 
Uniform weights • Mean vector x 

x is the reconstructed state in 
canonical units 


The mean 

variables 

0 = (*, % a, b) 

Mean sources s 
38 
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A planar eccentric orbit: 

• P 0 variances (1 km) 2 for position and (0.001 km/s) 2 for velocity 

• Ephemeris at true anomaly 0 = 0° after 1 period 

• N = 2000 particles 

• Period 18.1013 hrs. 

Initial Condition X 0 : 

[28000 km, 0 km, 0 km/s, 4.1331 km/s] 


Orbital 

Elements 

Numerical 

Value 

a 

42,000 km 

e 

0.2 

i 

0 rad 

Q 

0 rad 

0) 

0 rad 

M e 

0 

Period 

18.1013 hrs 
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Component PDF compression Legend 

Original 

and reconstruction: + compressed 


A 

FFT and WT to 
obtain 

coefficients i 

J 




A 

Compression 
rate of 50% 

J 


l 

( 

V 

Inverse FFT and 
WT to 
reconstruct 
PMF 

\ 

/ 


Component 

FFTD kl (o||c) 

WTD kl (o| | c) 

1 

0.1759 

0.0010 

2 

0.1003 

0.0188 

3 

0.0686 

0.1001 

4 

0.1121 

0.1568 





0.08 

0.06 

0.04 

0.02 

0 

-15 


-10 


FFT 008 

0.06 

0.04 

0.02 

0 


0.025 

0.02 

0.015 

0.01 

0.005 

0 

- 0.005 


3 
0.08 

0.06 

0.04 

0.02 

0 

-°^5 


0.03 


0.02 


0.01 


0 . 03 - 

0.02 

0.01 

0 

- 0.01 - 
-4 


0.025 r 
0.02 
0.015 
0.01 
0.005 
0 

- 0 . 005 .- 


WT 


- 0.01 

-4 


0.1 

0.08 

0.06 

0.04 

0.02 

0 


Component 1 










§§ 

m 

BSiiiiilllii 

TERMS # TERMS WT 

# TERMS FFT 

Cost of storing number of terms 


Coefficients c 

273 

260 

Wavelet Filters 

8 

N/A 

Abscissa PDF 

12 

12 

Mixing Matrix A 

16 

16 

Mean(s) = x 

4 

4 

ft 

80 

80 

a 

20 

20 

% 

80 

80 

b 

4 

4 


80 

80 

V a 

20 

I 20 


80 

80 

V b 

4 

4 


681 660 


Method 

# Terms 

Cost 

(Method/PF) 

PF Predictions 

8000 

1 

WT Compressions 

681 

0.0851 

FFT Compressions 

660 

0.0825 



■ PF 

■ WT 

■ FFT 


However : need to quantify accuracies based 
on the state particle reconstructions 
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TOTAL TERMS 








TERMS 

# TERMS WT 

# TERMS FFT 

Coefficients c 

273 

260 

Wavelet Filters 

8 

N/A 

Abscissa PDF 

12 

12 

Mixing Matrix A 

16 

16 

Mean(s) = x 

4 

4 

/i 

80 

80 

a 

20 

20 

% 

80 

80 

b 

4 

4 

TOTAL TERMS 

497 

476 


Cost of storing number of terms 


Method 

# Terms 

Cost 

(Method/PF) 

PF Predictions 

8000 

1 

WT Compressions 

497 

0.0621 

FFT Compressions 

476 

0.0595 


100 


Q_ 



PF WT FFT 


■ PF 

■ WT 

■ FFT 


However : need to quantify accuracies based 
on the state particle reconstructions 
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KLD distance is calculated at different compression rates: [0.1 0.3 0.5 0.7 0.9] 



Compression rates Compression rates 


• The WT coefficients has a strong bias between the approximate and detail coefficients 

• WT performs better at lower compression rates compared to the FFT 
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1. Demonstrated the use of the Particle Filter: 

a) A nonlinear estimator 

b) Capable of incorporating non-Gaussian uncertainties during predictions 


2. Demonstrated the use of the Independent Component Analysis (ICA) and Principal component 
analysis for dimensional reduction 

3. Demonstrated the use of Nonlinear Factor Analysis (NFA) followed by FastICA to achieve 
uncorrelated states for PDF compression 

4. The Wavelet Transform and the Fast-Fourier Transform demonstrated as effective methods for the 
compression and reconstruction of univariate PDFs. 


POTENTIAL FUTURE WORK 

1. Use other smooth functions to represent the PMF of the components versus WT and FFT for 
compression 

2. Determine the limitations of NFA with respect to the number of particles (apart from 
computational cost) 

3. Other nonlinear decorrelation methods could be implemented ex. The Nonlinear 
Independent Factor Analysis (NIFA) that uses Gaussian mixtures for sources 
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Thank you. 




Extended Kalman Filter (LaViola 2003 et al) 

• Nonlinear models assuming Gaussian uncertainties 

• Linearization about a current mean (computes Jacobian) 
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Unscented Kalman Filter (Julier and Uhlmann 2004 et al) 

• Uses a series of weighted sample points to approximate the mean and covariance 
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Based on Bayes Theory 


Posterior Density 
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\x/ 

Particle Filter 
3 particles 
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kurt(w T x) = E {(wM) 4 } - 3[E {{w T x) 2 }] 2 -E {(w 7 x) 4 } - 3 1 1 w 1 1 4 
J(w) = E {(w T x) 4 } - 3 1 1 w 1 1 4 +FQ | w 1 1 2 ) 

F is a penalty term due to the constraint on w. The online learning algorithm has the form 

Mt + 1) = w(t)± M (t)[x(t)(w(t) T x(t)f-3 1| w(t) || 2 w(0+/(|| w(t) 1 1 2 MO] 

^ is the learning rate sequence, f is the derivative of F/2 

The first two terms in the brackets are the gradient of kurt(w r x) and the third term is the 
gradient of F(|| w || 2 ). 

The fixed points w of the learning rule takes the expectation of the learning rule and equating 
the change in weight to zero 

E{x(w T xf}-3 1| w|| 2 w+ f(\\w\\ 2 )w = 0 
w = scalar x (E {x(w T xf } - 3 1 1 w \ \ 2 w) 

FastICA 


w{k) = E {x(w(k - l) r xf } - 3 w(k - 1) 

Divide w(k) by its norm (scalar x w) 

| w{k) T w(k- 1) |— > 1 , let k = k+1 


56 




The prior assumptions $; and # are given as follows: 

The variance is parameterized to v = log o 


x ~ N(f(s, 6), diag(e 2V " )) 
s ~ N(0,diag(e 2r ’)) 
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C = C q + C p =E[\nq(s,d\%j\ + E{-\n p(s, d,x | H)] 

C; : Y.E[^q(s,(N))]+^E[\nq(d j )] = £ 
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The uncorrelated sources s', are binned into a 2 7 histogram then transformed into a 
normalized discrete PMF p for compression. 


FAST-FOURIER TRANSFORM 

• The FFT calculates the coefficients of the discrete PMF as follows: 

A — t/-ix*-i) 

c(k ) = 2^p(J)e n where (k = 1, n) p(j) is the discrete probability value 

j = i normalized from the histogram 


WAVELET TRANSFORM 




Use a wavelet filter (Daubechies 2) to extract terms that represent the function p 
as a function of approximate and detail coefficients 
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Data 
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1. Demonstrated the use of the Particle Filter: 

1. A nonlinear estimator 

2. Capable of incorporating non-Gaussian uncertainties during predictions 


2. Demonstrated the use of the Independent Component Analysis (ICA) and Principal component 
analysis for dimensional reduction 

3. Demonstrated the use of Nonlinear Factor Analysis (NFA) followed by FastICA to achieve 
uncorrelated states for PDF compression 

4. The Wavelet Transform demonstrated better compression results over lower compression rates 
compared to the Fast-Fourier Transform 

5. The Kullback-Leibler divergence and the Kolomogorov-Smirnov tests were used as Information 
measures to quantify the reconstructions 

6. Ephemeris compression and reconstruction were able to achieve great results for compression 
rates 

7. The reduced number of terms can be potentially used for cataloging accurate non-Gaussian 

distributed ephemeris fi1 




1. Implementation of the Particle Filter using adaptive number of samples for numeric 
efficiency and redundancy elimination. 


2. Include a higher fidelity force model to incorporate additional relevant perturbations 
such as atmospheric drag, solar radiation pressure, third body perturbations (Sun and 
moon), J2 effects etc. 

3. Develop a systematic way of determining the number of hidden neurons to use as well 
as the number of iterations during the training phase (possibly optimize based on 
constraints) 

4. Other decorrelation methods that do not assume Gaussian initial conditions as the NFA 
does, could be implemented ex. The Nonlinear Independent Factor Analysis (NIFA) 
(nonlinear counterpart of the ICA) 

5. Optimize the value n in the determination of the number of bins 2 n for the histogram 
generator 

6. Use other smooth functions to represent the PDF of the components versus a PMF 
before compression 

7. Determine the limitations of NFA, NIFA with the number of particles (apart from 
computational cost) 
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